Exact solution of a model DNA-inversion genetic switch with orientational control 
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DNA inversion is an important mechanism by which bacteria and bacteriophage switch reversibly between 
phenotypic states. In such switches, the orientation of a short DNA element is flipped by a site-specific recom- 
binase enzyme. We propose a simple model for a DNA inversion switch in which recombinase production is 
dependent on the switch state (orientational control). Our model is inspired by the fim switch in Escherichia 
coli. We present an exact analytical solution of the chemical master equation for the model switch, as well 
as stochastic simulations. Orientational control causes the switch to deviate from Poissonian behaviour: the 
distribution of times in the on state shows a peak and successive flip times are correlated. 

PACS numbers: 87.18.Cf,87.16.Yc,82.39.-k 



Reversible and heritable stochastic switching between two 
different states of gene expression is a common phenomenon 
among bacteria and bacteriophage, known as phase variation. 
Phase variation is often linked to pathogenesis, and may help 
bacteria survive fluctuating environmental conditions (e.g., a 
host immune system) ^ . An important molecular mechanism 
for phase variation is site-specific DNA inversion 12]. Here, a 
short piece of DNA (the "invertible element") is excised from 
the genome and reinserted (strand-by strand) in the opposite 
orientation by a site- specific recombinase enzyme binding to 
sequences at the ends of the invertible element. Different 
states of gene expression correspond to the two orientations 
of the invertible element ("switch states"). A well-known ex- 
ample is the fim genetic regulatory system, which controls the 
production of type 1 fimbriae in Escherichia coli. These fim- 
briae are important in uropathogenesis |3||] . In the fim system, 
the FimE recombinase is produced more strongly when the 
switch is in the on state than in the off state y, |5|] . This phe- 
nomenon is known as orientational control. 

In this Letter, we present a simple and general stochastic 
model for a DNA inversion switch with orientational control. 
We solve this model analytically, allowing us to determine the 
range of stochastic switching behaviour possible for this type 
of switch. We find that non-Poissonian behaviour occurs, re- 
sulting in a peak in the probability distribution of time spent 
in the on state and correlations between successive flips. Such 
non-Poissonian behaviour could have important effects on the 
population dynamics of switching microbes in changing en- 
vironments. One key parameter (the concentration of the re- 
combinase not under orientational control) controls the degree 
to which our model is non-Poissonian; this parameter corre- 
sponds to the main point of environmental regulation for the 
fim switch. In contrast, bistable genetic switches such aspos- 
itive feedback loops |^ or mutually repressing genes III S] in 
general show only Poissonian behaviour (exponential waiting 
time distributions and uncorrelated flips). 

The Model — Our model DNA inversion switch, illustrated 
in Fig.[Tl contains three elements: the invertible DNA element 
and two types of recombinase enzyme (Ri and R2). The in- 
vertible element has two possible orientations (the "on" and 
"off" states). These correspond to alternative patterns of gene 



expression, leading to different phenotypic states; however, 
we model here only the core of the switch and not its down- 
stream effects. The switch can be flipped between its two ori- 
entations by either of the recombinases. The concentration of 
recombinase R2 is assumed to be fixed, while the production 
of Ri depends on the switch state: Ri is produced only in the 
on state. This feature of the model constitutes its orientational 
control and leads to its non-Poissonian behaviour. The model 
is represented by the following reaction scheme: 
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Here, 6'on and 5'off denote the on and off states of the switch. 
Reactions (fTal) describe the production and decay of recombi- 
nase Ri : Ri is removed from the system with rate constant 
ki (due to cell growth and division, which we do not model 
explicitly), and is produced at rate /c2 only when the switch is 
in the on state. Reactions (fTbl) describe switch flipping. This 
may be catalysed by Ri with rate constants k^^ (on to off) 
and kf^ (off to on). We shall mainly consider here the case 
kf^ = 0. Recombinase R2 can also catalyse switch flipping. 
The concentration of R2 (fixed in our model) is not explicitly 
included in the reaction scheme, but is implicit through a de- 
pendence of the rate constants k^^^^^ on the R2 concentration. 

We note that mean-field, macroscopic rate equations corre- 
sponding to the above reaction scheme yield only one steady- 
state solution corresponding to the average switch state and 
average concentration of Ri. The underlying deterministic 
structure of the model is thus not bistable. In this sense, our 
model is fundamentally different from the bistable reaction 
networks presented in |0, 13, S] . 

Our model is inspired by the fim genetic regulatory sys- 
tem [3]. In analogy with fim, Ri in our model represents FimE 
while R2 represents FimB. Environmental stimuli such as nu- 
trient conditions and temperature act on the fim switch largely 
through changes in the level of FimB : this would cor- 
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respond to variation of our key parameters /c' 
our model is highly simplified in comparison to fim, in that 
it neglects cooperative and competitive recombinase binding 
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FIG. 1: (color online) A) Schematic illustration of the model. B) 
A typical simulation trajectory. The solid line represents the time 
evolution of the number n of i?i -molecules, the shading denotes the 
switch position and r indicates the duration of a period of the switch. 
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FIG. 2: (color online) Probability distribution for the number n of 
Ri molecules, with k^^ — kf^ = /c4, for different values of k4. The 
symbols show simulation results and the solid lines are the theoret- 
ical predictions. For the case where k^ = 10^, the (dotted) line is 
a Poisson distribution with parameter (1 + kT /kf)~^kT /ki. The 
other parameters are: ki = 1, k^ = 100, kf = 0.001 and kf = 0. 
The insets show the typical simulation trajectories for each distribu- 
tion, where the time is expressed in units of the mean period (r). 



and the effects of other DNA binding proteins. Our objec- 
tive in this work is not to model the details of the fim system, 
as other authors have done but rather to address gen- 

eral questions about the behaviour of this type of switch. Our 
model is designed to be as simple as possible while retaining 
the key features of DNA inversion and orientational control; 
its simplicity allows us to obtain analytical results and to ex- 
plore a wide range of parameter space. 

We simulated the reaction scheme ([B using a continuous 
time Monte Carlo scheme A typical trajectory is shown 
in Fig. [B where we plot the number n of Ri molecules and 
the switch state as functions of time. When the switch is in the 
on state, n increases (on average) towards a plateau value of 
/C2//C1, while in the off state n decays towards zero. We now 
obtain an exact analytical solution for the case where kf^ = 0. 
This case is relevant to the fim switch, where FimE catalyses 
almost exclusively on to off switching. In the following, all 
our analytical results will correspond to kf^ = 0, while sim- 
ulation results will be presented also for kf^ > 0. Analytical 
results for kf^ > will be published elsewhere. 

Steady State. — We first consider the statistics of n in the 
steady state and compute the long time, joint probability Ps (n) 
that the switch is in state s and there are n molecules of Ri . 
The system of birth-death equations for Pon{n) and Poff(^) 
becomes in the steady state 

(n + l)kiPon{n + 1) + k2Pon{n - 1) + kfpoff{n) 

= {nki + ^2 + nfeg^ + fer)Pon(n) , (2a) 

(n + l)kiPof{{n + 1) + nA:§>on(n) + fefponln) 

= {nki + kf)poff{n) . (2b) 

In order to decouple the above set of equations, we solve ([2al) 
for poff, then insert the result into ([2bb to give a decoupled 



equation for pon- Introducing the generating function G5(z) = 
z^Ps{n) (where s = {on, off}), the decoupled equation 
for pon reduces to a second order differential equation for Gon- 
Defining then a new variable u = Uz = k2z/{ki -\- kf^) — 
^1^2/(^1 + /^3^)^, the latter equation reads: 

uGM + (a - u)G',,{u) - bGon{u) = , (3) 

where a = 1 + ^1 + + kf)/{ki + and b = 
1 + kf^/ ki . Expanding the solution as a regular power se- 
ries (i.e. Gon{u) = ^^cimU^), one finds that Gon(^) = 
ao a, ii), where iFi is a confluent hypergeometric 

function and ao is an integration constant which can be 
determined with the normalisation condition J^^Pon(^) + 
Poff{n) = 1 (or equivalently Gon(l) + Goff(l) = 1). This 
result can be rewritten in terms of the original variable n as 

Ponin) = aQ— — ^^j^iFiih^n.a^n.uo) , (4) 
n! (a)n 

where {a)n = (^{a + 1) . . . (a + n — 1). An expression for 
Poff can be derived by inserting Eq.(|4|) into Eq. ([2al) . In Fig|2l 
we compare the result for p{n) = Pon + Poff to simulations, 
obtaining perfect agreement. Figure [2] also illustrates the ef- 
fects of the different timescales for switch flipping and pro- 
duction/decay of Ri. For small k^, switch flipping is slow 
compared to the rate of change of n. In this case, when the 
switch is in the on state, the number of recombinase has time 
to reach a plateau before the switch flips off. The on and off 
switch states are then each associated with a different value 
of n and the distribution p{n) is bimodal. In contrast, when 
/c4 is large, the switch flips back and forth much more rapidly 
than recombinase production or removal. Then, the fraction of 
time spent in the on state is kf^/ {kf^ + k^^), and p{n) tends 
to the Poisson distribution expected for a birth-death process 
with birth rate k^kf/{kf + kf) and death rate ki. 
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Flipping time distributions. — To determine how orienta- 
tional control affects switch function, we compute flipping 
time distributions. The flipping time T can be defined in 
two different ways. In the first scenario, which we caU the 
Switch Change Ensemble (SCE), we define T as the time 
spent in a particular switch state — for example, F^^^{T) 
is the probability distribution for the time between the mo- 
ment the switch enters the on state and the moment it flips 
from the on to the off state. In the second scenario, which 
we call the Steady State Ensemble (SSE), we start observ- 
ing the cell at a random moment and measure the time in- 
terval between this moment and its next flip into the other 
state. F^^^{T) and F^^^^{T) may be relevant to the response 
of a population of switching cells to a sudden environmen- 
tal change. They also correspond to an experiment where 
one measures the time until the next flip, for cells sampled 
in the steady state [|12|] . To compute these distributions, we 
define Fs{T\no) as the probability that the system begins at 
t = in the s state with no recombinase and flips for the 
first time at T. Note that for kf = 0, the off to on flip- 
ping process does not depend on Ri and is governed by kf^; 
thus Foff(T) = kf exp{-kfT) is independent of no (for 
both the SCE and the SSE). However, the on to off flipping 
rate is no dependent, so that we average over the ensemble of 
initial states (characterised by the probability Won(^) of hav- 
ing n recombinase at the start of our measurement) to obtain 
the flip time distribution Fon(T) = Won(no)Fon(T|no). 
For the SCE, the initial condition is taken just after a flip, 
which implies that ^^^^^(no) = Poff(^o)/G'off(l). For the 
SSE, the initial condition is sampled in the steady state, yield- 
ing W^^^{no) = Pon(no)/Gon(l). To compute Fon(T), we 
first define the survival probability /ion(^, t) that, at time t, the 
switch is in the on state with n Ri molecules, without hav- 
ing flipped, given the initial condition /ion(^o, 0) = Won(^o)- 
The evolution equation for Hon is: 

dtKnin, t) = in^ l)kihon{n + 1, t) + k2hon{n - 1, t) 

- {nki + /C2 + n/cg^ + /cf )/ion(n, t) . (5) 

Defining a generating function hon{z^t) = ^^z'^ho^{n^t), 
it follows that Fon(t) = —dthon{ljt). Eq.(|5]) reduces to a 
partial differential equation for hon{z, t) which has the initial 
condition hon{z, 0) = Xlno ^on(^o)^^°- Its solution is 

ho^[z t) = e"^'^^4"+fc2(l-fciron))gfe2ron(2-feiron)(l-e-*/"°") 

X /ion(/ciron + e-^/""-(z - feir),0) , (6) 

where Ton = {ki + k^)~^ . Fq^{T) can then be computed 
for the different measurement scenarios using ho^{x^{)) = 
Gs{x)/Gs{l) with s = off for SCE and 5 = on for SSE. 

Peak in the distribution. — Our results, illustrated in FigJH 
show a striking effect of orientational control for this model 
switch: for the SCE, we can obtain a peak in the flipping 
time distribution. Such a peak has been postulated for the 
Jim switch where it might imply that the switch tends not 
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FIG. 3: (color online) Occurrence of a peak in the flipping time 
distribution Fon^^(T), for A;4 = 0.1. In the shaded region the 
peak is present, and it vanishes outside this region. The dashed 
line shows the same result for /c4 = 0.25. As k4 increases, the 
range of parameters for which there is a peak decreases. The unit 
of time is set by ki (i.e. ki — 1). The insets show examples of 
F!n^{T) (soUd lines) and F^^^{T) (dotted lines). In the left inset, 
for {k2,kT,k4) = (0.2,4,0.1), F!^^{T) is peaked, while in the 
right inset, for (5, 1, 0.1), Fon^{T) shows no peak. 

to leave the on state before it has had time to synthesise fim- 
briae Time spent in the on state may also influence 
recognition by the host immune system. This peak in F^^^ is 
a consequence of the feedback between the switch state and 
the level of Ri . Once in the on state, the rate of on to off flip- 
ping increases with time as Ri is produced. In contrast, for 
a Poissonian switch, the rate of flipping is constant in time. 
For a peak to occur, the slope of F^^^ at the origin must be 
positive, which implies 

k2 - {kTf/kT - ih + 2kT) (no) - k^ (ng) > 0, (7) 

where (...) denotes an average using the weight Won- The 
l.h.s. of d?]) can be evaluated numerically using the exact re- 
sult Q to compute (no) and (no). We can then determine the 
regions of parameter space where a peak exists, as shown in 
the shaded region of Fig.[3]for the SCE. Our results show that 
the presence of a peak is favoured by large values of /c2 (strong 
production of Ri in the on state), and suppressed by very 
large values of k^^ (strong -mediated on to off switching) 
or by very small values of k^^ (switching dominated by Ri- 
independent mechanism). Likewise, when /C4 is increased, the 
range of values over which the peak exists is decreased, since 

-independent Poissonian switching tends to dominate. For 
the SSE, on the other hand, we did not find any parameter 
values where inequality ^ is verified. The peak in the SCE 
appears because the number of recombinase, and hence the 
flipping probability, is typically low immediately after enter- 
ing the on state and increases significantly thereafter. In con- 
trast, in the SSE one typically starts a measurement when n, 
and hence the flipping probability, is already high. This tends 
to suppress the peak in the SSE flipping time distribution. 

Correlated flips. — Another potentially important effect of 
the feedback between the switch state and the production of 
recombinase Ri may be to cause correlations in the waiting 
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FIG. 4: (color online) Correlation function 

[i^i^j) ~ i^i)'^] I [('^^^) ~ ('^^)^] between switch periods Ti 
and Tj, for several values of the rate A;3^^. These are simulation 
results for kx = 1, k'f - 5, = 1, k't = A:f = 0.1. 

times between successive flips (for example, a particularly 
short time before a flip might lead to subsequent flips occur- 
ring in quick succession). Such correlations might allow a 
population of switching microbes to "remember" the history 
of past environmental changes. We define the switch period 
Ti as the time from when the switch enters the on state from 
the off state for the i\h time, until it enters the on state for 
the (i + l)th time [cf. Fig.[Tl). In Fig.lH we plot simulation 
results for the correlation function for switch periods and 
Tj , as a function of j — i. For kf = 0, weak correlation is 
observed between subsequent periods Ti and r^+i. Correla- 
tions are weak because when kf = 0, the off to on switching 
process does not depend on Ri and is an uncorrected Poisson 
process. When kf ^ 0, correlations are much stronger and 
extended correlated sequences of flips emerge. 

Discussion. — We have presented a generic model for a 
DNA inversion switch with orientational control. By solv- 
ing the model analytically in the case kf = (relevant to the 
Jim switch), and using stochastic simulations, we have shown 
that this type of switch can display markedly non-Poissonian 
behaviour, including a peaked flipping time distribution for 
intermediate values of k^^ and, for kf > 0, correlated se- 
quences of flips. Non-Poissonian behaviour has beenpos- 
tulated to be a consequence of orientational control U 0]. 
The model presented here allows us to analyse the origins 
and effects of this behaviour in detail, and provides analyt- 
ical results which can be used as a basis for more complex 
models lO . Suggested evolutionary roles for orientational 
control include rapid response to environmental change 
as well as peaked flipping time distributions fsl [^. This 



study raises interesting questions about the consequences of 
non-Poissonian switching for population dynamics in chang- 
ing environments. Several models have been proposed for the 
growth of populations of switching cells in stochasticall y an d 
periodically changing environments (see, for example uM)- 
These models assume Poissonian switch flipping. The an- 
alytical solutions presented here should make it possible to 
extend such models to the case of non-Poissonian flips. Non- 
Poissonian switch flipping opens up the possibility that lin- 
eages of cells may 'remember' (in a statistical sense) the his- 
tory of their recent phenotypic states. This is likely to have im- 
portant consequences for models which include selection ac- 
cording to the fitness of different phenotypic states in chang- 
ing environments. We speculate that bacteria which adopt 
a non-Poissonian flipping strategy may be able to maximise 
the evolutionary advantages of using some knowledge of the 
likely future behaviour of the environment, combined with the 
benefits of a stochastic strategy as an insurance against sudden 
and unpredictable environmental changes. These avenues will 
be the subject of future work. 
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